ON PARTIAL REGULARITY OF STEADY-STATE 
SOLUTIONS TO THE 6D NAVIER-STOKES 
EQUATIONS 



HONGJIE DONG AND ROBERT M. STRAIN 

Abstract. Consider steady-state weak solutions to the incom- 
pressible Navier-Stokes equations in six spatial dimensions. We 
prove that the 2D HausdorfF measure of the set of singular points 
is equal to zero. This problem was mentioned in 1988 by Struwe 
[24] ■ during his study of the five dimensional case. 



1. Introduction 

In this paper we consider the incompressible steady-state Navier- 
Stokes equations in six spatial dimensions with unit viscosity 

(1.1) uVu- Au + Vp = f, divM = 0, 

in a domain i7 C M^. We are interested in the partial regularity of 
weak solutions {u,p) to (11 .ip . 

Although the problem of the global regularity of solutions to the 
time-dependent Navier-Stokes equations in three and higher space di- 
mensions is still widely open, many authors have studied the partial 
regularity of weak solutions. In his pioneering work [191 |20l [22], Schef- 
fer established various partial regularity results for weak solutions to 
the 3D Navier-Stokes equations satisfying the so-called local energy 
inequality. In 3D, the notion of suitable weak solutions was first in- 
troduced in a celebrated paper [Ij by Caffarelli, Kohn and Nirenberg. 
They called a pair consisting of velocity u and pressure p a suitable 
weak solution if u has finite energy norm, p belongs to the Lebesgue 
space L5/4, u and p are weak solutions to the Navier-Stokes equations, 
and they satisfy a local energy inequality. After proving an e-regularity 
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criteria for local boundedness of solutions, they established partial reg- 
ularity of solutions and estimated the Hausdorff dimension of the sin- 
gular set. They proved that, for any suitable weak solution {u,p), there 
is an open subset where the velocity field u is regular and they showed 
that the ID Hausdorff measure of the complement of this subset is 
equal to zero. In [16], F. Lin gave a more direct and simplified proof 
of Caffarelli, Kohn and Nirenberg's result. A detailed treatment was 
later given by Ladyzhenskaya and Seregin in [TB]. Recently, some ex- 
tended results have been obtained in a number of papers; see Seregin 
[IB], Gustafson, Kang and Tsai [H], and Vasseur [2S], Kukavica [2], 
and the references therein. 

Much fewer results are available in the literature for the 4D and 
higher dimensional time-dependent Navier-Stokes equations, in which 
case the problem is more super-critical. In \21\, Scheffer showed that 
there exists a weak solution u in x R"*", which may not necessarily 
verify the local energy estimate, such that u is continuous outside a 
locally closed set of x R"*" whose 3D Hausdorff measure is finite. 
In a recent paper [2], the first author and D. Du proved that, for any 
local-in-time smooth solution to the 4D Navier-Stokes equations, the 
2D Hausdorff measure of the set of singular points at the first potential 
blow-up time is equal to zero. We remark that, in terms of the method, 
the dimension four in [2] is critical due to the following reason. To the 
best of our knowledge all the existing methods on partial regularity for 
the Navier-Stokes equations share the following prerequisite condition: 
in the energy inequality the nonlinear term should be controlled by 
the energy norm under the Sobolev imbedding theorem. Actually, four 
is the highest dimension in which we have such condition: L^^L^ 
-^oo-^f ^ fiv6 or higher dimensions this condition fails. 

This paper concerns the partial regularity of weak solutions u to the 
steady-state Navier-Stokes equations (11. ip . In the literature, the most 
relevant paper on the subject is [24j by Struwe, in which he proved 
the following e- regularity result, which implies that weak solutions are 
regular outside a singular set of zero ID Hausdorff measure. 

Theorem 1.1 (Struwe [21]). Let Vt be an open domain in R^ and 
f G Lq{Q) for some q > 5/2. There exists an absolute constant Eq > 
such that the following holds true. If u E H^{Q; R^) is a weak solution 
to (II. ip which satisfies a generalized energy inequality, and if for some 
xo E fl there is Rq > such that 



J\x—xo\<j- 

then u is Holder continuous in a neighborhood of Xq. 
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The proof of Theorem 11.11 rehes on some techniques for proving reg- 
ularity for elhptic systems (cf. Morrey [T7] and Giaquinta [12]) and 
estimates for the hnear Stokes' system due to Solonnikov [23J. A simi- 
lar approach was used before by Giaquinta and Modica in [llj to study 
the steady-state Navier-Stokes equations in dimensions d < 4. Be- 
cause time corresponds to two space dimensions, in some sense the 5D 
stationary problem is similar to the 3D non-stationary problem. There- 
fore, dimension five is the smallest dimension for which the steady-state 
Navier-Stokes equations are super critical. We also note that Theorem 
II. II was improved by K. Kang [13j, in which partial regularity up to the 
boundary was established for a smooth domain i7 C M^. The existence 
of regular solutions to the steady-state Navier-Stokes in high dimen- 
sions have also attracted substantial attention. We refer the reader to 
[HI |25l m El O [71 13] and the references therein. 

In [21] Struwe raised the following interesting question: does the 
analogous partial regularity result hold in dimension six or higher. It 
seems to us that some arguments in [21] do not work in six or higher 
dimensions. In this paper, we give a positive answer to Struwe's ques- 
tion in dimension six. To be more precise, we shall prove the following 
regularity result with a sufficiently small constant ^o- Let Q be an open 
set in M^, / G L6,ioc(^), and let m be a weak solution to (11.11) satisfying 
a local energy inequality (12. 2p . Then if for some xq G there exists an 
Rq > such that 



then u is Holder continuous in a neighborhood of Xq. In particular, it 
follows that the 2D Hausdorff measure of the set of singular points of 
the weak solution u is equal to zero. 

Related to Struwe's work, our proof also uses some techniques which 
appeared in the study of the regularity theory for elliptic systems; but 
our approach is quite different from that in [21]. In particular, we do 
not use any estimate for the linear Stokes' systems. Roughly speaking, 
there are three steps in our proof. In the first step, we essentially 
follow the argument in [2], which in turn used some ideas in [T3] and 
^16j . The novelties are in the second and third steps. In the second 
step, we choose suitable test functions in the energy inequality and then 
use an iteration method to establish a weak decay estimate of certain 
scaling invariant quantities. In the last step, we successively improve 
this decay estimate by a bootstrap argument, and finally we use the 
elliptic theory to get a good estimate of the L3/2 norm of Vm, which 
yields the Holder regularity thanks to the Morrey lemma. 
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It is, however, unclear to us whether similar results hold in seven or 
higher dimensions. In fact, we believe that six is the highest dimension 
to which our approach can be applied; see Remark 13.81 Therefore, in 
some sense, our results are critical in terms of the dimension. 

To conclude this introduction, we explain some notation used in the 
sequel: M*^ is the d- dimensional Euclidean space with a fixed orthonor- 
mal basis. A typical point in M'^ is denoted by x = {xi,X2, ■■■,Xd)- As 
usual the summation convention over repeated indices is enforced. And 
X ■ y = Xji/j = ^^=1 XjUj is the inner product for x,y E M'^. The stan- 
dard Lebesgue spaces are denoted by Lp {p> 1). Various constants are 
denoted by in general and the expression A^ = A^( - ■ ■ ) means that 
the given constant A^ depends only on the contents of the parentheses. 



For summable functions p,u = (ui) and r = (rjj), we use the follow- 
ing standard differential operators 



S{xo,r) = {x eR^\\x-xo\=r}, 5(r) = 5(0,r), S = S{1), 
B{xo,r) = {x eR^\\x-xo\ <r}, B{r) = B{0,r), B = B{1). 



We also denote the mean value of a summable function as follows 



Here |y4| as usual denotes the Lebesgue measure of the set A. 

Let Xq be a given point in f2 and r > a real number such that 
B{xo,r) C Q. It is known that in the sense of distributions one has 

Ap = -Dij [uiUj) + div / 

= -Dij {{ui - [ui].^,^^r)iuj - {uj]xo,r)) + dlv /, in B{xo, r). 

This will hold for a weak solution to ( 11. ip . Now let r]{x) be a smooth 
function on supported in the unit ball B{1), < rj < 1 and t] = 1 
on 5(2/3). We consider the decomposition 

(2.1) p = p^^^r + Ko,r, in B{xo, r). 

Above Pxo,r solves the following Poisson equation 

^Pxo,r = -Dij{{Ui - Nxo,r)(Mi - [Uj]xo,r)Vxo,r) + div(/r7^.o,r) , 



2. Setting and main results 



u^i = DiU, Vp = (p^i), Vn = (uij), 
divu = Ui^i, divr = {rijj), Au = divVu. 



These are all understood in the sense of distributions. 
We use the following notation for spheres and balls 
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where T];j.g^r{x) '= ri{{x — XQ)/r). Then h^^^^ is harmonic in B{xQ,2r/3). 

We will omit the indices of p and h whenever there is no possibility of 
confusion. The following notation will be used throughout the article: 

A{r) = A{r,xo) = — / {ul"^ dx, 



E{r) = E{r,xo) = — \Vu\^ dx, 



C{r) = C{r,Xo) = — I \uy dx, 

JB(xo,r) 



D{r) = D{r,Xo) = ^ [ \p - [h^^rf^^ dx, 

^ JB{xo,r) 

F{r) = F{r,Xo)= I \f\^ dx. 

JBixn.r) 



lB(xo,r) 

Notice these objects are invariant under the natural scaling for (11. ip : 

u{x) — )■ \u{\x), p{x) — )■ }?p{\x), f{x) — !■ \^ f{\x). 

We will use these quantities to study the regularity of 6D steady-state 
suitable weak solutions to (11. ip . 

We say that a pair of functions {u,p) is a suitable weak solution to 
(II. ip in n if n G Hl^^iVl) and p G L3/2,ioc(^) satisfy (ll.ip in the weak 
sense and additionally the generalized local energy inequality holds for 
any non- negative test function ip G C^{VL): 

(2.2) 2 / \Vu\^i)dx< I \u\^^il) + {\u\^ + 2p)u-V^ + f ■ui)dx. 
Jn Jq 

The existence of regular solutions to the Dirichlet problem of (II. ip in 
four dimensions have been obtained Gerhardt [8], in five dimensions 
by Struwe [25], and in five and six dimensions by Frehse and Ruzicka 
[H [6]; see also [TOl |5l [71 13] for other related results. We observe that 
the regular solution constructed in [6] satisfies (12. 2p . On the other 
hand, the uniqueness does not hold in general unless some smallness 
assumption is imposed on the data. 

By the Sobolev imbedding theorem, for any functions u G Hi^^ and 
r > 0, we have the following inequality 
(2.3) 

r / r \ 3/2 / /. \ 3/2 

/ \ufdx<N[ / \Vufdx] +Nr'^[ / \ufdx] . 

JBir) \JB{r) J \J B{r) J 

This implies that C{r) is well defined for a suitable weak solution. 
Next we state the main results of the article. 
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Theorem 2.1 (e-regularity criterion in terms of E). Let Q be an open 
set in M.^, f G Lg iod^), and suppose that the pair {u,p) is a suitable 
weak solution to (11.11) in Q. Then there is a positive number Eq sat- 
isfying the following property. Assume that for a point xq & Q the 
following inequality holds: 

(2.4) limsupE(r) < ^o- 

Then u is Holder continuous in a neighborhood of Xq . 

Theorem 2.2 (e-regularity criterion in terms of C, D and F). Let Vl 

be an open set in M^, / G L^^iodP), and suppose that the pair {u,p) is 
a suitable weak solution to (11.11) in Q. There is a positive number Eq 
satisfying the following property. Assume that for a point Xq & Q and 
for some po such that B{xq,Pq) C Q the inequality 

(2.5) Cipo) + D{po) + F{po)<Eo 

holds. Then u is Holder continuous in a neighborhood of Xq. 

Theorem 2.3 (Partial regularity). Let Q be an open set in M.^, f G 
Lg ioci^), and suppose that the pair {u,p) is a suitable weak solution 
to (II. ip in fl. Then the 2D Hausdorff measure of the set of singular 
points in f2 is equal to zero. 

These results are in the spirit of the work of CafFarelh, Kohn and 
Nirenberg |r|, where it was proved that for any suitable weak solution 
u to the 3D time- dependent Navier-Stokes equations the ID Hausdorff 
measure of the set of singular points is equal to zero. For 5D steady- 
state Navier-Stokes equations, results of this type were obtained by 
Struwe [21]. 

Remark 2.4. The assumption of the external force / in Theorems 12.11 - 
12.31 may be relaxed. It should be possible to only assume that / in 
Lp,\oc for some p G (3,6) or alternatively to assume that / is in certain 
Morrey spaces. However, we do not intend to find such a minimal 
assumption of / in this paper. 

3. The proof 
We shall prove the main theorems in three steps. 

3.1. Step 1. In the first step, we want to control A, C and D in a 
smaller ball by the their values in a larger ball under the assumption 
that E is sufficiently small. Here we follow the argument in [2J, which 
in turn used some ideas in [15] and [16]. These first few estimates do 
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not use the equation f ll.ip . Furthermore, in this section we often write 
C{r) = C{r, xq) and similarly for A, D, E and F. 

Lemma 3.1. Suppose 7 G (0, 1), p > are constants and B{xo,p) C 
Q. Then we have 

(3.1) C(7P) < N[^-^Ey\p) +^-^A'/\p)E'/\p) +^'Cip)], 

where N is a constant independent of'-f, p and Xq. 

Proof. Denote r = 7p. By using the Poincare inequality and Cauchy's 
inequality, we have 

/ \u\'^dx= / {lul"^ - [\u\%g^p) dx + / [\u\\g^pdx 

JBixa.r) JBixn.r) JBixn.r) 



' B(xo,r) J B{xo,r) J B{xo,r) 

<Np \Vu\\u\dx+ (-] / \u\'^dx. 

Jb{xo,p) VP/ Jb{xo,p) 



This is furthermore bounded by 

1/2/ /■ , . xl/2 

'B{xo,p) ' ^JBixo,p) 



<Np( [ |VMpc/x)^^V / \u\'^dx\ 

^JB(xn.o) ' ^JBixn.o) ' 



+ ( - ) / \u\^dx 

Pj Jb{xo,p) 



<Np'A'/\p)( [ \Vu\'dxy^\(-X ( [ 

^Jb{xo,p) ' \PJ ^Jb(xo,p) ^ 

Due to the Sobolev inequality (I2.3p . we obtain 
! \ufdx<N (f \Vu\^dx\ 

JB{xo,r) \_\JB{xo,r) J 

\Jb{xo,p) J \PJ Jb{xo,p) 

The conclusion of Lemma 13.11 follows immediately. □ 

Lemma 3.2. Suppose 7 G (0,1/4] and p > are constants, and 
B{xo,p) C Q. Then for any Xi G -B(xo,p/4) we have 

(3.2) D{jp,x,) < N[j'/'D{p) +j"'E'/\p) +j-'F'/\p)], 

where N is a constant independent of'j, p, Xq, and Xi. 
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Proof. Denote r = 7/9. Recall the decomposition of p introduced in 
fl2.1l) . By using the Calderon-Zygmund estimate and the Sobolev- 
Poincare inequality, one has 

(3.3) / \p.,A^)f'dx 

Jb{xi,t) 



J B{xi,r) ' JW' 

<n( [ \Vufdxy^\N( [ Iff^^dx 

^JB(xi,r) ' ^JB(xi,r) 



Similarly, 

\Pxo,pf^dx<N([ iVwpdx)^^^ 

B{xo,p) ' ^Jb{xo,p) ^ 



(3.4) 



+ iV( / l/r/^rfx)'^'. 

B(xo,p) ^ 



Since /i^cp is harmonic in B{xQ,2p/?))., any Sobolev norm of /i^o.p ^ 
smaller ball can be estimated by any Lp norm of h^^^p in a larger ball. 
Thus, using the Poincare inequality, one can obtain 



B{xi,r) 



\h -\h ] 1^/2 r/r 

\"'Xq,p l"'Xo,p\xi,r\ '^•^ 

<Nr^/^ [ \VK,,pf/^ dx < Nr^^/^ sup 

J B{xi,r) ' B{x-i,r) 

Further using the estimates for harmonic functions and the inclusion 

B{xi,r) C B{xi,p/3) C B{xo,2p/3), the above is 

/ \ 15/2 p 

(3.5) <n(-] / \K,,p - [h,,,pU,p\'/^ dx 

\PJ Jb{xi,p/3) 
15/2 

< 



- / \p-[hxo,pU,pf^' + \Px,,pf'dx. 

P) JB{xo,p) 



^{^o,p) 

Combining (13. 4p and (13. 5 p together yields 

3/2 



(3.6) 



/ \p-[hxo,pU,rf^dx<N( [ iVul^dxX 

JB{xi,r) ^Jb(xo,p) ' 

15/2 



+ iv(/ \ffl^d:^^'\N[''-\ I \p-{h.,^x,,pf''dx. 

^JB{xo,p) ^ \PJ JB{xo,p) 
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Since px^^^r + h^^^r = p = Pxo,p + hxo,p in B{xi,r), by Holder's inequality 

(3.7) / I [^a:o,p]a;i,r [^xi ,r]a;i,r | ^ dx 

JB{xi,r) 

= Nr^Wh 1 -\h ] {^/"^ = Nr'^Wv 1 - f» 1 



<N [ \Pxo,p\''^'+\Px„r\''^'dx. 
•JB{xi,r) 

From dM]), (I32D, (O and ([MD, we get 

/ \p- [hxurU,rf^ dx<N( [ iVwp dxY^^ 

J B{xi,r) ^Jb{xo,p) ' 

+ N(f \f\<^/^dxy^\N(-) ' [ \p-[hx„pU,pf^'dx. 

^JB{xo,p) ^ VP/ Jb{xo,p) 

Finally, by Holder's inequality, the lemma is proved. □ 

Note that the following estimates use the equation (11. ip . or more 
precisely they use the generalized local energy inequality (12. 2p . 

Lemma 3.3. Suppose 9 G (0,1/2] and p > are constants, and 
B{xo,p) C Q. Then we have 

A{9p) + Ei9p) < N9'^ [C2/3(p) + C(p) + C^'\p)D'"\p) + F{p)] . 

In particular, when 9 = 1/2 we have 

(3.8) Aip/2) + E{p/2) < N[C''^{p) + C{p) + C'/\p)D'/\p) + F{p)]. 
Here N is a positive constant independent of 9, p and xq. 

Proof. Let r = 9p. By Holder's inequality, 

A{r) < C2/3(r) < N9-^C^''^{p). 

To estimate E{r), in the energy inequality (12. 2p we choose a suitable 
smooth cut-off function ip = ipi E C^{B{xo, p)) such that 

< ^1 < 1 in B{xo, p), ^1 = 1 in B{xo, p/2) 

(3.9) iWil < Np-^, iV^^Ail < in5(xo,p). 
By using (12.20 and because u is divergence free, we get 



2 I \u\ dx 

P Jb{xo,p) 

-[ {\u\^ + 2\p-[h]x,J)\u\dx+ [ \u\\f\dx 

P Jb{xq,p) JB{xo,p) 
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Due to the Holder inequality and Young's inequality, one obtains 
/ \u\^dx < ( / \ufdxf\ I dxf'^ < p^C^'^p). 

Jb{xo,p) Jb{xo,p) Jb{xo,p) 

And 

\p - [h]xoJ\u\ dx 

B{xo,p) 

<(/ \p - [huj'/' dxY^\ I \u\UxY" 

J B(xo,p) J B{xo,p) 

< Np'D^I\p)C^'\p). 

Furthermore 

I \u\\f\dx<\ [ \u\^dx + p^ [ \f\^dx. 

JB{xo,p) P Jb{xo,p) Jb{xo,p) 

Then, collecting these estimates. Lemma 13.31 thus follows. □ 

As a conclusion of this subsection, we obtain 

Proposition 3.4. For any small Eq > 0, there exists ei = ei{eo) > 
small such that for any Xq € satisfying 

(3.10) limsupE(r) < ei, 

we have 

(3.11) A{po) + E{po) + C{po) + D{po) < Bo, 
provided that po is sufficiently small. 

Proof. For a given xq E Q satisfying (13.101) . choose pi > such that 
-B(xo,pi) C Q. Then for any p G (0,pi], by using (13. 8p and Young's 
inequality 

A(7p) + Ei^p) < N[C'/'{2^p) + C(27p) + /^(27p) + F(27p)]. 

This estimate, (13.11) and (13.21) . with 7 G (0, 1/8), together with Young's 
inequality again implies 

(3.12) A(7p) + E(7p) + C(7P) + D{-fp) 

< N i^'C'/'ip) + ^'l'D{p) + 7^C(p) + 7'A(p)] 
+ Ar7-5° [E{p) + E\p) + F(p)] + N-i^ 
< iV7' \A{p) + E{p) + C{p) + D{p)] + N^^ 

-rNr"' [E{p) + E\p) + F{p)\. 
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Since / G Le.ioc, by Holder's inequality, we have 

(3.13) F{p) < 11/11 

It is easy to see that for any £o > 0, there are sufficiently small real 
numbers 7 < 1/ a/STV and Si such that if (l3.1Up holds then for all small 
p we have 

+ iV7-50(E(p) + E\p) + F{p)) < £0/2. 
By using fl3.12p . we reach 

A{^p) + E{^p) + C{^p) + D{^p) < i [A{p)+E{p) + C{p)+D{p)] + '-^, 

which together with a standard iteration argument gives flS.lip for some 
Po > small enough. □ 

3.2. Step 2. In the second step, first we will estimate the values of A, 
E and D in a smaller ball by their values in a larger ball. Note that in 
this subsection all of the quantities implicitly depend upon the point 
Xi as A{r) = A{r,xi) unless it says so otherwise. 

Lemma 3.5. Fix constants p > 6 E (0,1/3] and B{xi,p) C fl. 
Then we have 

(3.14) A{ep) + E{ep) < Ne^A{p) 

+ Ne-WA{p) + E{p)f^ + D{p)) + Ne~^F{p), 
where N > is independent of p, 9 and Xi. 

Proof. We prove the lemma by using a suitably chosen test function in 
the generalized local energy inequality (12.20 . Let r = Op. We define 

i>2{x) = (r^ + \x - Xi\'^)~^, 

which clearly satisfies A-?/'2 = — 24r^(r^ + |x — Xip)~^ so that 

(3.15) AV^2 < in < -cr~^ in B{xi,p), 

for some constant c > independent of r. 

In the energy inequality (12. 2p we choose ip = ipiip2, where ipi is taken 
from (13.90 in the proof of Lemma [3.31 with the center xi in place of Xq. 
Then we have 

(3.16) - / \u\'^'4jiA'ilj2dx + 2 / |Vu|ViV'2C^a; 

Jb{xup) Jb{xi,p) 

< [ {|up(^/'2A^/'i + 2V^Ai- VV'2) 
+ (|u|2 + 2{p - [/i],,,p))u ■ (^iVV'2 + ^2VV'i) + / ■ uV'iV'2} dx. 
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After some straiglitforward computations, from (13. 9p and (I3.15p . it is 
easy to see tlie following properties: 

(i) For some constant c > 0, on B{xi,r) it holds that 

'ipi'ip2 = V'2 > cr~^, -ilJiAip2 = -Aip2 > cr"*^. 

(ii) In B{xi,p), we have 

\^iH < Nr-\ \^|JlV^p2\ + |^/'2V^i| < A^r-^ 

iV'sA^il + |V^i ■ VV^al < Np-^. 

These properties together with (I3.16p . the Young and Holder inequali- 
ties, yield 

(3.17) A{r) + E{r) < N[e^A{p) + e-\C{p) + D{p)) + Q-^'F^p)]. 
Owing to the Sobolev inequality (12. 3p . one easily gets 

(3.18) C{p)<N[A{p) + E{p)f\ 

Upon combining (13.17^ and (13.180 . the lemma is proved. □ 

Lemma 3.6. Suppose p > is constant and B{xi,p) C Q. Then we 
can find a 6i E (0, 1) small, where 6i does not depend upon p, such that 

(3.19) A{e,p) + E{e,p) + D^'\e,p) < \ [A{p) + E{p) + D'/'^p)] 

+ iV(^i) [A{p) + E{p) + D^'^{p)\ + N{e,) [F{p) + fV2(p)] , 
where N is a constant independent of p and Xi . 
Proof. Due to ([S^D and (I3A4D . for any 7, ^ G (0, 1/4], we have 

D^/\-iep) < N[-i^D^/\ep) + -i-^F^/^{ep) + -i~^E{ep)\ 

< N-f^-^D^/\p) + N-i-^F^/^{p) + N-i-^e^A{p) 

(3.20) + Nr^e-^[A{p) + E{p) + D^'\p)]"^ + iV7"2r6F(p), 
and from (I3.14p we have 

(3.21) AiriOp) + EiriOp) < N{-fefA{p) 

+ iV(7^)-3[A(p) + Eip) + D2/3(^)]3/2 ^ iv(^^)-6F(p). 

Now we choose and fix 9 sufficiently small and 7 = 6^^^ such that 
N[-f^e-^ + 7"'^' + (7^)2] < Ar^2/5 < ^/g^ 
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Upon adding f l3.20p and f l3.2ip . we obtain 

A{^dp) + E{^dp) + D^/^j9p) 

< \ [A{p) + D'l^p)] + N[A{p) + E{p) + D''\p)f' 

+ N[F{p) + F'l\p)l 

where N depends only on 9 and 7. After putting 9i = 7^, the lemma 
is proved. □ 

In the next proposition we will study the decay property of A, E, C 
and D as the radius p goes to zero. 

Proposition 3.7. There exists Sq > satisfying the following property . 
Suppose that for some Xq E Q and po G (0, 1) satisfying B{xo,po) C Q 
we have 

(3.22) C(po, xo) + D{po, xo) + F{po, xq) < Eq. 

Then we can find N > and ao G (0, 1) such that for any p G (0, po/8) 
and Xi G B{xq, po/8) , the following inequality will hold uniformly 

(3.23) A{p,xi) + E{p,xi) + C^/'^{p, xi) + D'^^^{p,xi) < Np"'^ , 
where N is a positive constant independent of p and Xi. 

Proof. Fix the constant 9i from Lemma 13.61 Due to (13. Sp . (13. 2p and 
(I3.22p . we may first choose e' > then = Eo{e') > sufficiently small 
such that, for any Xi G B{xo, po/8), 

A(po/4, Xo) + ^(po/4, Xo) + D^/^{po/8, x,) < ^, 

and 

(3.24) N{9,)V^ < 1/4, N{9,){eo + e'J') < e'/2. 
where N{9i) > is the same constant from (I3.19p . By using 

5(xi,po/8) C B(xo,po/4) C n, 

we then have 

ifiipo) := A{po/8, xi) + ^(po/8, Xi) + D^/\po/8, xi) < e'. 

By using (I3.24p and (13.190 with p = po/8 we obtain inductively that 

¥.(^^Po) = A{9lpo/8,x,) + Ei9'lpo/8,x,) + D^/%91po/8,x^) < e'. 

(Holding for A; = 1, 2, . . .). It then similarly follows from (13.240 and 
(I3A9|) that 

(3.25) ^{9',po) < Ividt'po) + N,{9t'pof. 
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Above, thanks to (13.131) . we have used the estimate 

F(^^Vo/8,Xi) + FV2(^^-ipo/8,xi) < iVi(||/|U,(B(.o,,o/2)): 



rvo)^. 



Now we use a standard iteration argument to obtain the Holder con- 
tinuity of 93. We have to be a bit careful however because we do not 
make the standard assumption that ip{p) should be a non-decreasing 
function. We iterate (13.251) to obtain 



(3.26) ¥.(^fpo)<Q) ^(Po) + iVipg|^Q) 



< - 



^(Po) + 



2iVi 2' 



In the last inequality, without loss of generahty we have used that 
61 e (0, 1/2]. Since p G (0,po/32) we can find k such that O'lf < 4p < 
Then 

A{p,xi) + E{p,x,) + D^/%p,x,) < iV(0l)((/.(^^Vo)+i^'/'(^^Vo,Xl)), 

where we used Lemma 13.21 to estimate the third term on the left-hand 
side. By (I3.26P and (13.131) . the above is further bounded by 

k r 



Ni9^ 



m 2 

— TTPo 



+ N{ei)p' < Np' 



In this last line = A^(6'i, v9(po), A^i, po) and = ^jj^ 
yields (I3.23P for the terms A, E and D. The inequality for C{p,Xi 
follows from (13.181) . 



> 0. This 



□ 



3.3. Step 3 — Proofs of Theorems 12.1112.31 In the final step, we 
are going to use a bootstrap argument to successively improve the 
decay estimate (I3.23p . However, as we will show below, the bootstrap 
argument itself only gives the decay of E{p) no more than p^, i.e. one 
can obtain an estimate like 



\Vu\ 



V£ > 0, 



^dx < N{e)p^-', 

B{xi,p) 

for any p sufficiently small. Unfortunately, this decay estimate is not 
enough for the Holder regularity of u since the dimension is six (so that 
we need the exponent 4 + e instead of 4 — e according to the Morrey 
lemma). Then to fill in this gap we will use the elliptic theory. 

First we prove Theorem 12.21 We begin with the bootstrap argument. 
We will choose an increasing sequence of real numbers {afcjfcLi ^ (c^O; 2) 
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such that for any small 5 > we can find an integer m = m{6) with 
the property that > 2 — 6. 

For a fixed 6 > and m = m{S), under the condition (12. 5p . we claim 
that the following estimates hold uniformly for all p > sufficiently 
small and Xi G B{xo, Pq/8) over the range of {ak}^=Q- 

(3.27) Aip,xi) + E{p,xi)<Np''\ 

C{p, xi) < Np^'"'/^ D{p, xi) < A^p3"fe/2_ 

We prove this via iteration. The k = case for f l3.27p with ao was 
proven in f l3.23p . 

We first estimate A{p,xi) and E{p,xi). Let p = 6p where 6 = p^, 
p = p^~^ and p G (0, 1) to be determined. We use Lemma ESI and then 
f l3.27p (for ak) to obtain 

A{p) + E{p) < N (^p2M+afe(l-M) + p|afc(l-M)-3/. ^ ^4(1-,.)-6m ^ _ 

Choose p = -I "'^ , then fl3.27p is proven for Aip) + Eip) with the 
exponent of 

ak+i := min |2/i + ak{l - p), ^«fc(l - - 3p,4(l - p) - 6;u| 

3 12 

= - /i) - 3/i = — afc G (afc, 2). 

2 lU + Ofc 

Then the estimate in fl3.27p (with ak+i) for C{p, xi) follows from ( I3.18p . 
To prove the estimate in ( I3.27P (after level k) for D{p,xi) we will use 
Lemma [3. 2[ From 03.20 we obtain 

The estimate used here for E{p) follows from (13.130 . Now for any r 
small, we take the supremum on both sides with respect to p G (0,r) 
and get 

sup D{-fp,xi) < A^7^/^ sup D{p,xi) + N'j-^r^°"^+^/'^ + N'j-^r^. 

pe(o,r] pe(0,r] 

Since 9/2 > 3 > fa^+i, by using a well-known iteration argument, 
similar to f l3.25p (or see e.g., \12\ Chap. 3, Lemma 2.1]), we obtain the 
estimate in f l3.27p (with ak+i) for D{p). Then we have shown how to 
build the increasing sequence of {ak} for which fl3.27p holds. Moreover, 

10 , , 10 , 
10 + Qfe 10 + ao 
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which imphes that — )■ 2 as A; oo. Note that by the above proof, 
the constant in fl3.27p may go to infinity as A; — ?■ oo; thus we truncate 
at level m < oo. 

In particular, f l3.27p with k = m gives for any small 6 = 6{m) > 
that 

6-<5 

5 



(3.28) / \u\^dx<Np 

Jb{xi,p) 

(3.29) [ \u\^ + \p-[h]^.,,pf^dx<Np^ 

We obtained these estimates via the bootstrap argument, next we will 
use the elliptic theory to improve them. 

Now we fix a 5 G (0, 1/10) and rewrite (11. ip (in the sense of distri- 
butions) into 

Aui = DjiuiUj) + Dip - fi. 

Finally, we use the classical elliptic theory to complete the proof. 
Thanks to (I3.28p . there exists pi G (p/2,p) such that 

(3.30) f \u\^dx < Np^-^. 

J S{xi,pi) 

Let V be the unique solution to the Laplace equation 

Avi = inS(xi,pi), 

with the boundary condition Vi = Ui on 5'(xi,pi). It follows from the 
standard estimates for harmonic functions, Holder's inequality, and 
(ICTl) that 



(3.31) sup \Vv\<Npf \v\dx < Np-'^-^/^ 

B{xi,pi/2) Js{xupi) 



Denote w = u — v E H^{B{xi, pi)). Then w satisfies the Poisson 
equation 

Aw^ = Dj{uiUj) + Di{p- [/i]xi,p) - fi in 5(xi,pi). 

with zero boundary condition on 5'(xi,pi). By the classical Lp esti- 
mates for the Poisson equation, we have 

iivt/;iu3/,(B(.„p,)) < ^iii^nL3/.(B(.„po) 

+ N\\p- [hU,ph^^^^six„p,)) + Npi II/IIl3/,(b(x,,pi)) • 
This together with the assumption on / and fl3.29p gives 
(3.32) ||V«;|U3/,(5(.„p,)) < Np'-' + Np' < Np'-'. 



PARTIAL REGULARITY OF 6D STEADY SOLUTIONS TO NSE 



17 



Since |Vm| < \Vw\ + |Vf |, we combine fl3.3ip and f l3.32p to obtain, for 
any r G (0,p/4), that 



[ I Vm I '^/^ dx < Np^-^^/^ + Nr 

JB(xi,r) 

Upon taking r = p^/'^'~^l^ jA (with p small), we get 
(3.33) / \Vu\^/^ dx <Nr^, 

JB(xi,r) 



6^-3/2-3(5/4 



where 



5/4-5/8 



Since ( 13.33P holds for arbitrary xi G B{xo, po/8) and all r small, by 
the Morrey lemma (see for instance fl2\ Theorem 1.1 on p. 64 of Ch. 
Ill]), u is Holder continuous in a neighborhood of Xq. This completes 
the proof of Theorem 12. 2[ 

Theorem 12. II then follows from Theorem 12. 2l by applying Proposition 
13. 4[ Finally, Theorem 12.31 is deduced from Theorem 12.11 by using the 
standard argument in the geometric measure theory, which is explained 
for example in [1], or alternatively in |12j . 

Remark 3.8. Finally we remark that by using the same method we can 
get an alternative proof of Theorem 1 1.1 1 for the 5D steady-state Navier- 
Stokes equations if we assume that / G Lsjoc- However, it seems to 
us that six is the highest dimension to which our approach (or any 
existing approach) applies. In fact, by the Sobolev imbedding theorem, 
if^(M^) L3(]R^). So the nonlinear term in the energy inequality can 
be controlled by the energy norm when d = 6 but not higher. 
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